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ABSTRACT 

We consider the primordial gravitational wave (GW) background in a class of spatially- 
flat inflationary cosmological models with cold dark matter (CDM), a cosmological 
constant, and a broken-scale-invariant (BSI) steplike primordial (initial) spectrum of 
adiabatic perturbations produced in an exactly solvable inflationary model where the 
inflaton potential has a rapid change of its first derivative at some point. In contrast 
to inflationary models with a scale-free initial spectrum, these models may have a GW 
power spectrum whose amplitude (though not its shape) is arbitrary for fixed ampli- 
tude and shape of the adiabatic perturbations power spectrum. In the presence of a 
positive cosmological constant, the models investigated here possess the striking prop- 
erty that a significant part of the large-angle CMB temperature anisotropy observed 
in the COBE experiment is due to primordial GW. Confronting them with existing 
observational data on CMB angular temperature fluctuations, galaxy clustering and 
peculiar velocities of galaxies, we find that for the best parameter values f?A ~ 0.7 
and h w 0.7, the GW contribution to the CMB anisotropy can be as large as that of 
the scalar fluctuations. 

Key words: cosmology: theory - early Universe - cosmic background radiation - 
large-scale structure of Universe. 



1 INTRODUCTION 



There seems to be increasing observational evidence in sup- 
port of the inflationary paradigm (see the review in Linde 
1990; Kolb & Turner 1990; Lyth & Riotto 1998), which of- 
fers an elegant solution to some of the outstanding problems 
of standard Big Bang cosmology. During inflation, primor- 
dial quantum fluctuations (Hawking 1982, Starobinsky 1982, 
Guth & Pi 1982) of some scalar field(s) (inflaton(s)) are pro- 
duced, which eventually form galaxies, clusters of galaxies 
and the large-scale structure of the Universe through gravi- 
tational instability. 

However, the simplest inflationary CDM model in flat 
space (fi = 1) and with an approximately flat (ns ~ 1) ini- 
tial spectrum of adiabatic perturbations does not agree with 
observational data because it has too much power on small 
scales when normalized to the COBE data at large scales. 
One of the simplest and most elegant ways to reach agree- 
ment with observations is to add a positive cosmological 
constant to dustlike CDM. On one hand, this improves pure 



CDM models a great deal; on the other hand, the cosmolog- 
ical constant itself is viable in the presence of cold dark mat- 
ter only (as was emphasized e.g. in Kofman & Starobinsky 
1985). Now, the ACDM model is perhaps the most promis- 
ing CDM variant (Bagla, Padmanabhan & Narlikar 1996; 
Ostriker & Steinhardt 1995). A recent argument in favour 

of n < i, (n = i - n A ) (and n ~ (0.2 - 0.4)) follows 

from the evolution of rich galaxy clusters (Bahcall, Fan & 
Cen 1997; Fan, Bahcall & Cen 1997; Eke et al., 1998). The 
most recent data on SNe la events at large redshifts also 
strongly favour < 1 (Perlmutter et al. 1998; Garnavich et 
al. 1998; Riess et al. 1998), while the latest CMB constraints 
seem inconsistent with these low Q values in case Qa = 
(Lineweaver 1998). Finally, galaxy abundance at large red- 
shifts also disfavours open models (Peacock 1998), leaving 
the ACDM model with ns ~ 1 as the most successful model 
for structure formation. 

On the other hand, the observational fact that the 
smoothed slope of the initial spectrum ns cannot be sig- 
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nificantly different from 1 does not exclude the possibility 
of local strong deviations from scale invariance, i.e., steps 
and/or spikes in the initial spectrum. Of course, such a be- 
haviour is not a typical one, and may be expected to oc- 
cur only in exceptional points in Fourier space (if at all). 
At present, one scale k w 0.05/iMpc _1 is a candidate for 
this role, since there exists an evidence for a peculiar be- 
haviour (in the form of a peak or a sharp change in ns) in the 
Fourier power spectrum of rich Abell-ACO clusters of galax- 
ies around this point (Einasto et al. 1997a, 1997b, 1997c; 
Retzlaff et al. 1998). This evidence is, certainly, inconclusive, 
especially since there is no peak at k « 0.05ftMpc _1 in the 
power spectra of both APM clusters (which are generally less 
rich than Abell-ACO clusters) and APM galaxies (Tadros, 
Efstathiou & Dalton 1998), though a break in ns may ex- 
ist in the latter spectrum too (Gaztanaga & Baugh 1998). 
Note that the attractive possibility to explain this feature by 
Sakharov (acoustic) oscillations fails (Atrio-Barandela et al. 
1997, Eisenstein et al. 1997). Hence, if confirmed in the fu- 
ture, this feature should originate in the initial perturbation 
spectrum itself. 

Therefore, there exist serious reasons to continue the in- 
vestigation of BSI CDM models. Then, it is natural to fit the 
arbitrary wavenumber of the singularity appearing in these 
models to k = 0.05ftMpc _1 . From the theoretical point of 
view, inflationary models with BSI spectra having such a 
characteristic scale naturally appear as a result of a phase 
transition approximately 60 e-folds before the end of infla- 
tion (Kofman & Linde 1986; Adams, Ross & Sarkar 1997; 
see also Starobinsky 1998). However, it is important that in- 
flation cannot produce very sharp features in the primordial 
spectrum. In particular, an exact ^-function-like step (which 
would result from a naive application of the slow-roll expres- 
sions to the case in which the first derivative of an inflaton 
potential V(ip) sharply changes near some point >p = <po) is 
not permitted; what appears instead follows from the exact 
solution found in Starobinsky (1992). 

Summarizing, if we believe both in the present cosmolog- 
ical constant and in the cluster data exhibiting a preferred 
scale in the perturbation spectrum, we have to consider the 
ACDM model with a BSI initial spectrum. Remarkably, a 
detailed comparison with the bulk of all existing observa- 
tional data shows that the ACDM model plus the above 
mentioned BSI spectrum gains much from the cosmologi- 
cal constant too (Lesgourgues, Polarski & Starobinsky 1998, 
hereafter Paper I). Not only does the inclusion of A > en- 
larges the allowed region of cosmological parameters fi and 
Ho, but it also permits the new possibility of an inverted 
step (i.e. more power on small scales) in the initial spec- 
trum. Also, this model avoids the main problem found in a 
model of double inflation (Lesgourgues & Polarski 1997) for 
which the Doppler peak turns out to be low even for those 
values of the parameters with an acceptable matter power 
spectrum P(k). However, in Paper I, the free parameters of 
the model were chosen in such a way that the primordial 
GW background generated during inflation in addition to 
adiabatic (scalar) perturbations was too small to contribute 
significantly to the observed large-angle CMB temperature 
anisotropy. Now we want to investigate if it is possible to ob- 
tain a large GW contribution to the temperature anisotropy 
AT/T in this model. 

Historically, the primordial GW background generated 



from quantum vacuum metric fluctuations during inflation 
was the first observational prediction of inflation: its spec- 
trum was first calculated in Starobinsky (1979) even before 
the first viable models of inflation were constructed, and the 
CMB temperature anisotropy for the multipoles I = 2, 3 pro- 
duced by this background was found by Rubakov, Sazhin & 
Veryaskin (1982); see Lyth & Riotto (1998) for references on 
numerous subsequent papers. The GW and adiabatic con- 
tributions to AT/T cannot be separated if only the latter 
quantity is measured, so the quantity of practical interest 

is actually \J 1 + — 1, i-e., the relative excess of the ob- 
served rms value of AT/T over the rms value calculated un- 
der the assumption that there is no GW contribution at all. 
However, measurement of the CMB polarization provides a 
unique opportunity to distinguish both these contributions 
and to prove the existence of a primordial GW background 
directly. In the case of the simplest inflationary models with 
scale-free spectra, Cf /Cf can be very small, e.g., in case 
of the "new" inflation or the R + R 2 inflation (where R 
is the Ricci scalar). On the other hand, for chaotic infla- 
tionary models with a power- law inflaton potential V(tp), 
Cio /Cfo is comparable with unity, though still rather small, 
e.g., Cf /Cfo « 0.2 for the quartic potential. Another exam- 
ple is provided by power-law inflation with an exponential 
potential V(ifi) leading to a tilted CDM model. However, 
such models did not prove to be successful in explaining the 
height of the first accoustic (Doppler) peak for ns < 0.9, 
while ns > 0.9 results in (Cfo/Cfo) < 0.5. Note that infla- 
tionary models having non-negligible Ci /Cf may be also 
characterized by the condition that the total variation of 
the inflaton field during the period corresponding to present 
scales in the range (10 - 10 4 )/i _1 Mpc is not negligible com- 
pared to the Planck mass (Lyth 1997). 

The obstacle for having a sufficiently large CJo/Cfo in 
the case of a scale-free initial perturbation spectrum lies in 
the observational fact that the bandpower 1(1 + 1)C; grows 
almost by an order of magnitude when I changes from 10 
to 200 - 300 (see Table 1 below). This is even larger than 
what is expected in the standard ns = 1 CDM model with- 
out GW and a cosmological constant, while the presence of a 
significant GW contribution to AT/T on large angles should 
manifest itself in a decrease of the height of the first Doppler 
peak relative to the normalization at I = 10. In other words, 
a high first Doppler peak is an argument against signifi- 
cant primordial GW background in case of scale-free initial 
perturbation spectra. Note that even using old CMB data 
from the Saskatoon 94 and South Pole 94 experiments re- 
ferring to I = (70 ± 20), it was already possible to reach the 
conclusion that {C[ /Cf ) < 0.7 with 97.5% probability in 
the cases of chaotic and power-law inflation with h — 0.5 
and Q = 1 (Markevich & Starobinsky 1996). Now, with the 
whole set of data presented in Table 1, this upper limit be- 
comes much lower. The existence of a positive cosmological 
constant only slightly relaxes this argument, so that we still 
get (Cj /Cfo) < 0.15 in this case (see Sec. 3.2 below, the 
case p = 1). 

Thus, at present, the only possibility to have a signifi- 
cant primordial GW background is to require some kind of 
significantly non-scale-invariant initial density power spec- 
trum. The first attempt in this direction was performed by 
Lukash & Mikheeva (1996, 1998) who considered the case of 
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the inflaton potential V(tp) = Vq + m £ with m 2 not small 
as compared to Hq = StvGVo/3. However, this model has 
a "blue" initial spectrum with ns > 1 at small scales, and 
faces serious problems regarding its excess of power at scales 
smaller than 8h~ Mpc if ns > 1.3. On the other hand, if the 
parameters of the model are taken in such a way that the 
asymptotic slope satisfies ns < 1.3, then (CTo/Cfo) < 0.5. 
So, it appears that in order to get (CJo/Cfo) = 1 or more, 
one has to take a BSI spectrum with a significantly sharper 
feature in it. This gives one more reason to investigate the 
model considered in Paper I with respect to the possible 
existence of a large primordial GW background. 

Therefore, as in paper I, we suppose that the inflaton 
potential derivative, V'(ip), has a rapid change from A+ to 
A—, when tp decreases, in a neighbourhood of tpo. This gen- 
erates a steplike spectrum of primordial adiabatic fluctua- 
tions, with a specific substructure. In a vicinity of the step, 
the amplitude of large-scale and small-scale plateaus is given 
by: 



81 HE 



8nG 



(1) 



where $ is the (peculiar) gravitational potential at the mat- 
ter dominated stage. So, the amplitude of the step is given 
by p = A-/A+, and we call ko the characteristic scale of the 
step (for details, see paper I). In this model, it is still pos- 
sible to fix freely the amount of primordial GW for given p 
and normalization, since their initial spectrum near k = ko 
is given, for each polarization state, by the standard expres- 
sion: 



k 3 (h mn h mn )(k) = IQ-kGHI 



(2) 



which does not depend on A± in contrast to k 3 & 2 (k) (here 
m, n = 1, 2, 3). Without the inclusion of a cosmological con- 
stant, we would be forced to consider the case p > 1 only, in 
order to increase power on large scales. Since A > already 
produces a desired excess of large-scale power, we are now 
free to consider both cases p > 1 and p < 1. 

However, in order to confront the model with observa- 
tional data, we need the expression for the primordial spec- 
tra far from the point k — ko, up to \y\ — |ln(fc/fco)| ~ 5. 
Since we adopt the natural assumption that any deviation 
from the slow-roll regime requires some special origin (e.g., 
some kind of phase transition during inflation) and, thus, 
that it should be an exceptional phenomenon, we assume 
that both slow-roll conditions given above are valid every- 
where far from the point k — ko (before the end of inflation, 
of course). Then standard expressions for the initial spec- 
tra of adiabatic perturbations and GW are obtained in this 
region. So, the full initial spectra for these perturbations 
follow from (|l|) and (@), respectively, by the substitution: 

2 _ 8irGV k 
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where the index k means that the quantity is taken at the 
moment of the first Hubble radius crossing (k — aH) during 
the inflationary stage. 

Note that since the second derivatives of V(tp) are not 
fixed by Eq.(l), we may freely take ns(k) far from the break 
point k — ko. However, \ns — 1| should be small in this 
region due to the two slow-roll conditions mentioned above. 
In paper I, we assumed that the smooth part of \ns — II is 



so small that it can be neglected at all for | In(fc/feo)| < 5, 
so the upper and lower plateaus of the scalar spectrum may 
be considered as flat ones. This assumption is self-consistent 
in the case of a negligible GW background. In our case, the 
situation is more complicated. In the slow-roll regime we 
have the well-known relation: 
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where the latter expression implicitly defines tpk as a func- 
tion of k. The coefficient 5 appearing here is approximate, 
its exact value depends on both h and SI (see Bunn, Liddle 
& White 1996 for a fitting expression at I = 14, which can 
be easily transformed to I = 10). Further, small corrections 
proportional to nr and (ns — 1) should be also taken into 
account if one wants to make this coefficient even more ac- 
curate (see Lidsey et al. (1997) for a review). Note that if we 
neglect the existence of a radiation dominated stage in the 
past (which corresponds to the approximation t eq <C tree), 
it would be equal to 5.3 for fl = 1, ns = 1, tit = (as was 
given in Polarski & Starobinsky 1995). 

Since we are interested in the case when Cfo /Cf is not 
small, we cannot assume that nr ~ 0. So, the amplitude 
of the GW background determines the slope of its power 
spectrum. Moreover, the values of nr on the right and on 
the left sides of the break point are different: 



nr{k > ko) — nr(k < ko) 



A 2 



8nGV 2 



(•») 



Strictly speaking, this equation refers to the smooth part of 
nr defined by differentiation of H 2 only; in addition, there 
is some substructure near k = ko- However, for our choice 
of ko, the only quantity of interest is nr(k) for k = (0.01 — 
0.5) feo- 

Therefore, we have to specify more accurately our initial 
spectra (or, equivalently, V(ip)). As a result, the BSI CDM 
model which we consider in this paper is actually different 
from the one with negligible GW studied in Paper I (though 
they have the same structure at k w fco). We shall consider 
the two most representative cases. 

1) ns fn 1 for both positive and negative large values 
of ln(fc/fco). In the slow-roll regime, this corresponds to an 
inflaton potential approximately proportional to (tp — l p±)~ 2 
for p > ipo and tp < po respectively, where ip± are two 
different constants. As discussed above, the choice of the 
initial scalar spectrum with flat plateaus at large and small 
fc's is very suitable for an explanation of both the small-scale 
observational data and the COBE results. Then, however, 
we cannot assume nr(k) to be constant: it must satisfy the 
equation (the "consistency relation" for single-field slow-roll 
inflation) which takes the following form for ns — 1: 

dxiT 



din k 



n T 



(6) 



(see, e.g., Eq.(14) in Polarski & Starobinsky 1995 and the 
more general discussion in the review by Lidsey et al. 1997). 
As expected, this case leads to the largest possible values of 
Cio/ Cio- 

2) ns ~ 1 for positive large values of ln(fc/fco), but 
nr — ns — 1 = const < for negative large ln(fc/fco). This 
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corresponds to the exponential form of V(tp) above tpo ■ Here 
we may assume constant slopes far from the break point 
k = ko- 



2 CONFRONTATION WITH OBSERVATIONS 

We want now to compare our model with observations. We 
use experimental data which constrain the matter power 
spectrum P(k) on one hand and the radiation power spec- 
trum on the other hand. As we have already noted earlier 
(Lesgourgues & Polarski 1997) when dealing with a model of 
double inflation, another model with broken scale-invariant 
primordial spectrum, the constraints on both types of fluc- 
tuations are essentially complementary, so that successfull 
confrontation of one spectrum of fluctuations does not pre- 
clude bad results for the other spectrum. One has to remem- 
ber that the point of this analysis is to investigate how large 
the amplitude of the produced GW can be while leaving the 
model still in good agreement with observations. Hence, we 
must adopt a strategy that enables us to answer this ques- 
tion with sufficient accuracy while leaving the possibility to 
refine the analysis later on if required by new observational 
evidence. 

The primordial scalar spectra under investigation have 
four free parameters: an overall normalization factor Qio, p, 
ko, and Cf /Cf . Given these four parameters, and some 
cosmological parameters, the scalar and (running) tensor 
tilts are completely set, for each possibility ns = 1 or 
ns — 1 + tit- We choose the step location as well as the 
overall normalization fixed: ko = 0.016 h Mpc -1 , which cor- 
responds to a bump in the present matter power spectrum 
P(k) at k = 0.05 frMpc~\ while we take Q 10 = 18uK (Ben- 
nett et al. 1996). The parameter p is found by requiring that 



locity at R — 50 h 1 Mpc, with Gaussian smoothing radius 
R s = 12/i" 1 Mpc: 

290 < V 50 < 460 km s _1 (1 - a confidence level) (7) 

in the absence of cosmic variance, and: 

245 < Vfso < 505 km s _1 (8) 

when cosmic variance is taken into account. As we will see, 
the lower bounds implied by these results yield stringent 
constraints on our models. 

• the STROMLO-APM redshift survey which gives a 
count-in-cells analysis of large scale clustering. Like other 
redshift surveys, it probes the matter perturbations well into 
the linear regime. Results are given in Loveday et al. (1992) 
for cells of nine different sizes. We compare these data points 
with the power spectrum (normalized to as = 1.00, and 
convolved with the nine corresponding window functions) 
through a \ 2 analysis. Since we can vary 3 parameters, there 
arc 6 degrees of freedom. 

• the power spectrum of rich Abell galaxy clusters, con- 
sisting of 36 points taken from Einasto et al. (1997a). Since 
a few points on the largest and smallest scales suffer from 
large uncertainties, we perform a \ 2 analysis with only 30 
points, corresponding to 0.023 < k < 0.232 AMpc -1 . As al- 
ready said in the introduction, this spectrum exhibits a clear 
feature with a bump at k = 0.05 h Mpc -1 and fco is chosen so 
that the matter power spectrum also exhibits a feature at the 
right scale. Hence, we are testing here a theory with 4 free 
parameters, and the \ 2 distribution has 26 degrees of free- 
dom. Strictly speaking, the Einasto et al. (1997a) data are 
not independent. Taking into account their correlation will 
decrease the effective number of degrees of freedom, how- 
ever this does not noticeably change our conclusions. We do 
not assume any specific value for the biasing factor for these 



a 8 = 0.60 fr ' 56 (White, Efstathiou & Frenk 1993; see also dato: for each set of parameters (h, Ov, C 10 /C 10 ), we cal- 
Viana & Liddle 1998), where as is the variance of the to- 
tal mass fluctuation in a sphere of radius 8/i _1 Mpc. Then, 
all models are simultaneously normalized to COBE (large 
scales, scalar plus tensor components) and to as (small 
scales, scalar component only). Note that p is still a function 
of the remaining three free parameters: the two cosmological 
parameters h and J1a on one hand, the inflationary param- 
eter CJo/Cfo on the other hand. In this way, the latter is 
singled out as the only remaining inflationary free parame- 
ter. 

We then compute numerically, using the fast Boltzmann 
code CMBFAST by Seljak & Zaldarriaga (1996), the matter 
power spectrum P(k) and the CMB power spectrum C(ff) 
for various values of the parameters h, 0,a , Cfo/ Cfo . We will 
find it interesting to state the results in two-dimensional cuts 
of the parameter space for several given values of h. This is 
particularly appropriate if we expect other observations to 
yield some refined a priori knowledge of the value of h. 

In order to constrain the matter power spectrum, we 

use: 



culatc and adopt the biasing factor yielding the smallest \' 
In that sense, this test probes only the shape of the power 
spectrum. 

In order to constrain the radiation power spectrum, we 



• the CMB anisotropy data on large and small angular 
scales found in various experiments. We use the bandpower 
estimates AT; ± a given for an experiment characterized by 
a window function Wi. More precisely, we have: 



AT, 



a bs(T) 

Vim ' 



(9) 



where a bs(T), the observed rms temperature fluctuation 
(which differs from the theoretical value by the inclusion of 
the window function Wi characterizing the experiment), is 
given by: 



1 



1 = 2 



c,w. 



(10) 



• peculiar velocities taken from the MARK III cata- 
log (Willick et al. 1997) and the POTENT reconstruction 
(Bertschinger & Dekel 1989; Kolatt & Dekel 1997) of the 
velocity field. Peculiar velocities have the advantage that 
they probe all mass (and not just galaxies), but they have 
rather large uncertainties. We use here the rms bulk ve- 



and I(Wi) is defined as: 

W) = D' + 5)i(iTi)- 

1=2 

The bandpower estimate assumes a flat spectrum 1(1 + 
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1)C oc C2 being a good approximation over the width of 
the window function Wi centered around some "effective" 
multipolc number l e defined through l e = I (IWi) / 1 (Wi) . 
We compare the theoretical CMB power spectrum with 19 
experimental points given in Table 1, again through a \ 2 
analysis, with 16 degrees of freedom^]. The Saskatoon data 
provide an important constraint on any candidate model, as 
they give a clue to the correct height of the first Doppler 
(or acoustic) peak. However these are given in Netterfield 
et al. (1997) with a calibration error of ± 14%, which ap- 
plies equally to all five points. This remaining uncertainty 
is treated as in Lineweaver & Barbosa (1998a, 1998b): for 
each set of free parameters (h, 0,a, Cfo/Cf ), we make a 
preliminary \ 2 analysis in order to find the best calibration. 
Then, the Saskatoon data are treated on equal footing with 
other CMB data, and a second \ 2 is computed with the 19 
points. 

We chose to perform three separate \ 2 analysis for the 
Stromlo-APM, the CMB and the rich Abell cluster data. 
In this way, we see clearly what is implied by each set of 
data separately. We further avoid to put on the same footing 
observations of a very different kind. 



3 RESULTS 

The result of all the tests are given in Figure [l], in both case 
ns = 1 ans ns — 1 + bt, for a few two-dimensional cuts of 
the parameter space corresponding to h = 0.5,0.6,0.7,0.8. 
For peculiar velocities, we plot the curves corresponding to 
the lower limit V50 = 245 km s _1 (the limit without cosmic 
variance, V50 = 290 km s _1 , is also given for completeness). 
For the three \ 2 tests, we plot contours inside the preferred 
regions, in which % 2 i s smaller or equal to the number of 
degrees of freedom. Let us comment these results step by 
step. Since the results for the ns = 1 and ns — 1 + jit are 
quite similar, we will discuss them simultaneously. 

3.1 Peculiar velocities and STROMLO-APM 

First, let us not consider the rich Abell galaxy clusters dis- 
tribution. In this case, the matter power spectrum is con- 
strained, on one hand, by bulk velocities, and on the other 
hand, by the STROMLO-APM redshift survey. The bulk 
velocity test probes the amplitude of P(k) at intermediate 
scales, corresponding to the step in the primordial spec- 
trum. This amplitude can vary significantly, even if erg has 
been fixed on smaller scales, through the matter transfer 
function and the primordial spectrum (indeed, a change in 
Cio/Cfo induce a change in p). Hence, this test depends 
on all three parameters. One can see on Figure [l] that for 
(ns = 1, h < 0.6), and for (ns = 1 + «t, h < 0.8), this 

* After this paper was submitted, the results of two new impor- 
tant CMB experiments have been published. They confirm the 
data presented in Table 1, and, therefore, strenghten our con- 
clusions. The first experiment, QMAP (de Oliveira-Costa et al. 
1998) confirms the Saskatoon normalization at I = 90, and the 
shape of the Saskatoon data up to I = 150. The second one, the 
NCP at OVRO (Leitch et al. 1998) confirms the second CAT 1 
point at I K, 600. 



constraint cuts the preferred region in the Cfo/Cfo - fl A 
parameter space below the value Cfo/Cfo = 1. 

On the other hand, the STROMLO-APM redshift sur- 
vey probes the shape of P(k) only at scales smaller than 
the step, i.e. those scales given by the high-fc plateau of the 
primordial spectrum. Hence, it only depends on the mat- 
ter transfer function, not on Cfo/Cfo, and on Figure |l| it 
appears as a constraint on Qa only. 

These two constraints define a wide preferred region in 
parameter space, and it is crucial to include CMB anisotropy 
measurements in order to obtain good predictions. 

3.2 CMB anisotropics 

The CMB x 2 test selects preferred elliptic regions in the 
Cio /Cf — planes, which truly admit an extension up to 
Cio/Cfo = 1 an d even more. This is due to the 'inverted' 
step in the primordial spectrum, which compensates for the 
loss of power in small-scale anisotropy usually implied by 
a high tensor contribution. The GW contribution even im- 
proves the model: in all cases plotted in Figure [j], the best 
values are obtained for 0.4 < Cfo/Cf < 0.8. 

For all h considered here, intersections with the previ- 
ously preferred regions are found, and enhanced in yellow (or 
gray) on Figure El As expected from what was just said, the 
major part of the preferred regions correspond to 'inverted' 
steps 0.45 < p < 1. Those with p = 1 (Harrison-Zeldovich) 
and 1 < p < 1.1 ('usual' step) correspond to the lower left 
corners of the preferred regions (not shown in Figure 2), 
with Cfo/Cfo < 0.15. 

3.3 Clusters distribution 

Let us now consider the rich Abell galaxy cluster distribu- 
tion data which constitutes the main motivations for the 
characteristic scale appearing in our model, and which de- 
termines its location through the value of ko- The elliptic 
preferred regions (see Figure |l|) are narrow when compared 
with STROMLO-APM regions, because data error bars are 
smaller, and sensitive to larger scales (where p, and there- 
fore Cio/Cfo values are crucial). The best values are ob- 
tained in the absence of gravitational waves, but values of 
Cfo/Cfo — 1 are still acceptable (except for ns = 1 + nr, 
h < 0.6). There is no intersection with the CMB preferred 
region when h < 0.5. The intersection of all preferred re- 
gions is coloured in bright yellow (or dark gray) on Figure 
H, showing that: 

• for h > 0.7, larger preferred region in the Cf /Cf — 
SIa plane are obtained, with 0.68 < Q,a < 0.77 and < 
Cfo/Cfo < 1. 

• for h = 0.6, a rather narrow region survives, for which 
the GW contribution is substantial, Q,a ~ 0.63 and 0.2 < 
C10/C10 < 1. 

• for h = 0.5, a marginal intersection between the CMB 
and cluster favoured regions is found for 57a ~ 0.53. This 
possibility is still acceptable, because we are rather severe 
in the definition of the preferred windows. 

In the previous discussion, we made almost no distinc- 
tion between ns — 1 and ns — 1 + nr, since the results are 
quite similar in both case. This shows that our predictions 
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Table 1. Bandpower estimates used for the CMB \ 2 test. 
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do not depend very much on the assumptions on large scale 
tilts, provided that the consistency relation is satisfied. In 
particular, the results would still hold for an arbitrary con- 
stant scalar tilt in the range 1 + nr(ko) < ns < 1- The 
only difference is that the preferred parameter windows and 
the lowest \ 2 regions are systematically obtained for smaller 
Cf /Cf when ns = 1+nr, while the first case, with ns ~ 1, 
admits slightly larger values of Cf /Cf . 



4 CONCLUSION 

The generation of a GW cosmological background is an es- 
sential prediction of all inflationary models. Amplitude and 
statistics of this background can be computed from first 
principles, as also their contribution to the CMB temper- 
ature (and polarization) anisotropies. However, only part of 
the inflationary models have a primordial GW background 
sufficiently large in order to be actually measured, using 
either CMB measurements, or else gravitational-wave an- 
tennas for their direct detection. An accurate measurement 
of the CMB angular temperature anisotropy and the CMB 
polarization (especially on angular scales > 1°) may result 
in the remarkable discovery of a primordial GW background 
on cosmological scales which, of course, would be of great 
importance by itself, and would also provide a very strong 
argument for the inflationary scenario of the early Universe. 
The already measured height of the first Doppler peak for 
multipoles I = 200 — 300, in spite of all the indeterminacies, 
is so high that it precludes a significant GW contribution 
to the multipole dispersion Cio (where Z = 10 is the char- 
acteristic multipole for the COBE data) if the initial power 
spectrum of adiabatic perturbations is scale-free. Shifting 
to a BSI initial spectrum helps to avoid this obstacle and to 
keep the possibility of finding the GW background. 

In this paper, as in Paper I, we have considered the BSI 
CDM model with a positive cosmological constant and we 



have confronted it with recent observational data on CMB 
temperature anisotropies, large-scale galaxy-galaxy correla- 
tions, peculiar velocities of galaxies, and spatial correlations 
of rich Abell clusters. Both new ingredients of this model, 
as compared to the standard CDM model with a flat initial 
perturbation spectrum, namely the cosmological constant 
and the particular form of the BSI initial spectrum, have 
been already introduced earlier to account for hitherto un- 
explained data, without reference to a primordial GW back- 
ground. Now we have shown that this model, with a slightly 
different choice of its parameters as compared to the one 
considered in Paper I, admits a large GW background, too. 
A positive cosmological constant is essential for this since 
the initial power spectrum admits an inverted step, p < 1 
(i.e., more power on small vs. large scales), without which 
the quantity (Ct/Cs)io characterizing the relative GW con- 
tribution to Cm would be small. 

We have considered two subclasses of this model, differ- 
ing by the behaviour of the scalar (density) power spectrum 
in its approximately scale-free part far below the break- 
point at A; = ko- In the first case, the spectral slope is 
ns ~ 1 in this region (while we may not assume neither 
rir ~ 0, nor neglect its scale dependence). In the second 
case, riT = ns — 1 = const < 0. For each class, we have 
investigated how large this GW background can be while 
still leaving the model in good agreement with observations. 
The difference between the two cases is small, as expected. 
We find that the best models are located in essentially the 
same window of the cosmological parameters h and Qa: 
h ~ 0.7, S1a ^ 0.7 while the GW contribution to the tem- 
perature fluctuations on large angular scales, measured here 
c T 

by the parameter -^2-, can be as large as that of the scalar 

fluctuations. Note that this is still of course many orders 
of magnitude below a value allowing their possible direct 
detection by ground-based gravitational-wave antennas. 
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Figure 1. Results of all the tests, in both case ng = 1 ans ng = 1 + nj, for a few two-dimensional cuts of the parameters space 
corresponding to h = 0.5,0.6,0.7,0.8. The preferred regions for each test are limited by dotted curves in blue (bulk velocities), red 
(STROMLO-APM), black (cluster distribution), and green (CMB anisotropy). Inside the former three regions, we plot the contours 
associated with integer values of \ 2 - 
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